The dependence of transverse and longitudinal resolutions on incident Gaussian beam 
widths in the illumination part of optical scanning microscopy 
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We studied both theoretically and experimentally the intensity distribution of a Gaussian laser 
beam when it was focussed by an objective lens with its numerical-aperture (NA) up to 0.95. 
Approximate formulae for full widths at half maximum (FWHM) of the intensity distribution at 
focus were derived for very large and very small initial beam waists with respect to the entrance 
pupil radius of the objective lens. In experiments the energy flux through a 0.5 micron pinhole was 
measured for various pinhole positions. We found that the FWHM's at focus in the transverse and 
the longitudinal directions do not increase much from the ultimate FWHM's until the input beam 
waist is reduced below the half of the entrance pupil radius. In addition, we observed significance of 
the spatial distribution of the input beam against a true Gaussian beam profile in the case of small 
initial beam waist. For high NA with resulting focal beam waists comparable to or smaller than 
the wavelength of the laser, the interaction between the electric field and the conducting surface 
of the pinhole caused the transverse FWHM to be measured slightly smaller than FWHM of the 
unperturbed intensity distribution convoluted with the pinhole opening. 

PACS numbers: 



I. INTRODUCTION 

The spatial resolution in optical scanning microscopy 
is critically dependent on both the beam spot size near 
the focus of a scanning objective lens and how the focal 
spot is imaged back onto an imaging plane. In order to 
achieve the ultimate resolution the beam spot size at fo- 
cus should be minimized for a given illumination source. 
Otherwise, the spatial resolution is degraded and it can- 
not be recovered however well one handles the imaging 
of the focal spot. For proper accessing the minimal focal 
beam spot, one should be able to calculate and measure 
the beam spot size accurately. 

In many experiments using an objective lens we usually 
assume that the incident beam is a plane wave apertured 
by the entrance pupil of the objective lens. However, the 
light source in the optical scanning microscopy is often a 
Gaussian laser beam, not an ideal plane wave. One can 
expand the Gaussian beam and let the central part of it, 
simulating a plane wave, incident on the objective lens. 
A practical question is then how large the beam should 
be expanded with respect to the entrance pupil size of 
the objective lens in order to obtain a spatial resolution 
comparable to that with the ideal plane wave input. 

To answer this question, we need to know the near- 
focal plane intensity distribution of a Gaussian laser 
beam with an initial beam width wq when focused by an 
objective lens with an entrance pupil diameter D 
The intensity distribution, in general, can be calculated 
by the electromagnetic diffraction theory of Richards and 
Wolf HH|| . This theory is based on the vectorial equiv- 
alent of the Kirchhoff-Frcsncl integral in the Debye ap- 
proximation 0,0- 
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The intensity distribution in the region of focus have 
been measured in several experiments by using a knife- 
edge jHH [Hill and a tapered fiber However, 
a systematic investigation of the near-focus intensity dis- 
tribution in the non-paraxial regime as a function of the 
input Gaussian beam width wq has not been reported. 

In the present work, we re-examine the diffraction the- 
ory of Richards and Wolf for input Gaussian beams. In 
two limiting cases of very small and very large input beam 
widths, we derive approximate formula for the full width 
at half maximum (FWHM) of the intensity distribution 
at focus in the longitudinal and transverse directions. We 
then confirm the validity of theoretical predictions in ac- 
tual experiments employing objective lenses with numer- 
ical apertures of 0.4, 0.75 and 0.95 for various Gaussian 
input beam widths. 

This paper is organized as follows. In Sec. [H] we 
first theoretically examine transverse and longitudinal 
FWHM's near the focal plane for an arbitrary input beam 
waist wq and then derive approximate formulae for limit- 
ing cases, Wq <C R and Wq ^> R. Experiment is described 
in Sec. II I II and results and discussion are presented in 
Sees. HVl andlVl We summarize the work in Sec. IVII 



II. THEORY 

Suppose a Gaussian beam with a waist wo is incident 
on an objective lens with a high NA and an entrance 
pupil radius of R. We can think of three different regimes, 
namely, (i) wq <C R, (ii) wq ~ R, and (iii)wo 3> R- We 
first consider a general theory which can address all three 
regimes and then discuss regimes (i) and (iii) as limiting 
cases of the general theory. 
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FIG. 1: Coordinate system for the calculation of the intensity 
distribution in the region of focus. 



A. Field distribution near the focal region in 
general cases 

We use the electromagnetic diffraction theory of 
Richards and Wolf 0, 0, 13 f° r the numerical calcula- 
tion of the intensity distribution of the focused beam. 
For the integral, we choose our Cartesian coordinate sys- 
tem in the following way (see Fig. 1). The origin is lo- 
cated at the focus, z axis coincides with the optic axis 
of the optical system under consideration, pointing in 
the beam propagation direction and x axis points in the 
polarization direction of the incident field eo. A time- 
independent part e(r) of the analytical solution of the 
Helmholtz's equation for the electric field at a point P(r) 
in the image space of our optical system is given by 1 



o-x = erj(#) •v/cosi? [cost? + sin 2 tp(l — cost?)], 

a y — eo(t?) v / cost?(cost? — 1) cos ip sin ip, 

a z = — eof^Vcos •& sin t? cos <p , (2) 

where eo(t?) is the amplitude of the incident electric field 
eo- Similar expressions hold for the components of the 
magnetic field strength factor b. On substitution of 
Eq.(J2| into Eq.QJ with s = (sin i? cos <p, sin t? sin ip, cos t?), 
we obtain the following expressions for the Cartesian 
components of e. 

i 

e x (r) = - -kf(I Q + I 2 cos 2(f)), 

1 

e y( r ) = ~-^kfI 2 sm2(j), 

e z {r) = -ik fix cos 0, (3) 

where 

J o (r,0) = / e (i?) %/cost? sin 0(1 + cost?) 
Jo 

x Jo(kr sin i? sine?) exp(ikr cost? cose?) til?, 

h{r,6) = / e (t?)v / cost?sin 2 i? 
J o 

x Ji(kr sin t? sin ff) exp(i£x cost? cos ff) d-d, 
h{r,0) = / e (t?)v / cos t? sin t?(l - cost?) 

xJ 2 {kr sin t? sin f?) exp(iA;r cost? cos ff) d(t4,) 

where a is a semi-aperture angle satisfying Q = 2ir(l — 
cos a) and its Sine value is the numerical aperture 
(NA=sina). 

For a well-collimated Gaussian beam with a beam 
waist wo and an amplitude Aq, erj(i?) can be written as 



e ( r )=-^T J J a( " Sx s ' Sy ^ exp{ik[<i>(s x , s y )+s-r}} ds x ds t 

where s = (s x , s y , s z ) is a unit vector pointing in the 
direction of a ray, $(s x ,s y ) represents aberration in the 
optical system, f2 is the solid angle subtended by the 
exit pupil of the objective lens from the origin, the focus, 
and a, called an electric strength factor, is the electric 
field incident on the exit pupil after passing through the 
lens. Similarly, the magnetic field h(r) can be written 
in the same way in terms of a different strength factor 
b(= s x a). Eq.QJ is valid only if kf 1, where / is the 
focal length. 

We introduce spherical polar coordinates (/, t?, <p) for 
the point Q on the exit pupil and (r, 9, <f>) for the ob- 
servation point P in the image space. The Cartesian 
components of the strength vector a can then be written 
as 



e (t?) = A,exp[-(/ sint?/ Wo ) 2 ]- (5) 

under the Abbe's sine condition [l5| . 

The quantity to be measured in our experiment to be 
presented below is the power transmitted by a small aper- 
ture near the focal plane. This quantity is nothing but 
the time-averaged z-component of the Poynting vector, 
which is given by 

S z (r) = ^(|/ | 2 -|/ 2 | 2 ), (6) 
where c denotes the speed of light in vacuum. 

B. Large beam waist limit, wo 3> R 

Since wq 3> R, we can approximate the incident Gaus- 
sian beam as a plane wave and use the results in the pre- 
vious section with a substitution eo(t?) = Aq (constant) 
in Eq. @J. 



3 



1. Transverse spot size fAiFWHiij 

The field distribution in the focal plane of the objective 
lens can be written as 



h{r,6 
h(r,6 

h(r,e 



tt/2) 
tt/2) 
tt/2) 



A / Vcos$sintf(l + costf) 
Jo 

x Jo(fcr sin ■&) d-d, 

A / Vcos -d sin 2 i9 
j o 



x Ji(fcrsinz?) di?, 



Ao / V cos sin — cos i?) 
Jo 

x J 2 (fcr sin z?) dd, 



(7) 



In general, I ^> Ii, I 2 and thus the transverse spot 
size at focus is mostly determined by Iq integral. Fur- 
ther approximation is then obtained by noting that the 
functional factor (1 + cos-#)/2 is approximately equal to 
y/cosd, which can be easily verified by Taylor series ex- 
pansion of these two. This approximation is reasonably 
good even when # ~ 1. For example, the difference be- 
tween these two fuctional factors is 4.8% for i? = 1. Under 
this approximation, Eq. (JJJ becomes 



Jo 



Ji(fcr sin a) 



2A / cos -d sin -d J (kr sin -d) d-d cx 

/ fcrsma 



which is of the same form as the Fraunhofer diffraction 
by a circular aperture. Although the paraxial assump- 
tion sin a <gc 1 is used in the Fraunhofer diffraction the- 
ory, our approximate result, Eq. ©, is still applicable 
to non-paraxial cases with a up to the order of unity. 
This finding is new and has not been recognized. The 
transverse spot size is then obtained from Eq. JSJ) as 



Ax 



2 x 1.6163 



FWHM 



0.5145 



A 

NA 



(9) 



Figure [5] shows the difference between Axfwhm approx- 
imated by Eq. © and the exact one by Eqs. (@J and JSJ- 
The approximation is excellent in that the difference is 
as small as 2.8% even when NA=1, the largest possible 
NA value. 




0.0 0.2 0.4 0.6 0.8 1.0 
NA 

FIG. 2: Error in Aifwhm approximated by Eq. @ with 
respect to the exact one by Eqs. Q and © as a function of 
NA. 



Under the same approximation as above, 
r a 

Iq w 2^4 / cos z9 sin $exp(«/cz cos 
Jo 

q exp(iq)dq 



{kz) j f, z co „ , , 



oc sin a 



smx 



i tan 



a x cos x — sin x 



(11) 



(8) where x = fczsin 2 (a/2). 



For a up to unity, the contri- 
bution from the second term in |io| 2 is negligibly small, 
proportional to tan 4 ~ < 0.089, and thus |/o| 2 is approxi- 
mately given by the Sine function squared, which is again 
the same as the Fraunhofer diffraction result except that 
x is proportional to a 2 not to sin 2 (a/2) in the usual 
Fraunhofer diffraction. For an arbitrary a, Azfwhm is 
obtained from Eq. (jTT|) as 



rj(a)X 7/(arcs 

AIFWHM = . . 2 a = a ■ 2 l\ 

4sur I 4snT (± 



?7(arcsinNA)A 
9 (iarcsinNA) 



(12) 



where the slowly varying function r)(a) is plotted in Fig. 
13 For a up to unity, we can approximate r/(a) ~ 77(0) ~ 
1.772, by which our error is only 1.7% for a = 1 and 5.4% 
for a — 1.25, which corresponds to NA=0.95. Under this 
approximation, 



2. Longitudinal spot size (Azfwhm) 

The field distribution in the z-axis near the focus is 
jiven by 



I o (r = z,0 = O) = A Q / \/costfsintf(l + cosi?) 



x exp(ikz cos -d) d-d, 
h(r,6 = 0) = h{r,6 = 0) = 0, 



Azfwhm 



1.772A 



1.772A 



(13) 



4 sin f 4sin (iarcsinNA) 
which reduces to the usual Fraunhofer diffraction result 



Azfwhm =s 1.772-^ ~ 1.772 



NA^ 



(14) 



(10) under the paraxial condition, a«l. 
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FIG. 3: Numerical factor 17(a) in Eq. fESI . 



C. Small beam waist limit, wo <C _R 



Although the numerical aperture of the lens is assumed 
to be large, only the central portion of the objective lens 
is utilized by the incident Gaussian beam when wq <C R- 
One can define an effective numerical aperture NA e g- as 
NA e fj = wq/ f <C 1, and thus the paraxial approximation 
can be effectively applied. One is allowed to use Gaussian 
optics to calculate the beam size in the focal region. Par- 
ticulary, when the incident beam has a minimum waist 
at the entrance pupil of the objective lens, the Gaussian 
optics provides a simple formula for the field distribution 
in the region of focus. 



1. Transverse spot size (AiFWHiij 



the focal plane can be written as 
I (r,6 = ir/2) = 



eo(^)vcos^sini?(l + cosd) 
x,J (kr sin??) d$, 



eo {'&) V cos?? sin 2 D 
x,Ji(kr sinz?) d/Q, 



eo 0?) V cos 1? sin i?(l — cost?) 



x J 2 (kr sin tf) dtf, (17) 

where eo($) is given by Eq. JSJ. Since eo(#) is significant 
only when sini? < wq/ f 1, the integrands above count 
only when and thus we can rewrite the above as 



J « 2 / eo(i?)i?Jo(fcn?)d0, 
Jo 

/•a 

Ji pb / e (??)tf 2 Ji(fcn?)dtf, 



h ~ ^ e (tf)i? 3 J 2 (fcn?)dtf, 



(18) 



Since Jx/Jo - (W/) 2 « 1 and J 2 /J - (W/) 4 « 1, 
the field distribution is mostly determined by I Q . We can 
further simply the Iq integral as 

I oc ex-p[-(ftf/w f}-dJo(krd)dtf 





rfa/wo 

x I exp(— x 2 )x J( 





( krw 



exp(-x 2 )xJ (px) dx = exp[-(p/2) 2 ](19) 



where p = krwo/f, from which we obtain an 1/e width 
of the field distribution as 2//fcw , which is nothing but 



•j' in Eq. (JUJ. 



The Gaussian beam waist wq' in the region of focus is 
given by 



w 



TTWq 



(15) 



where wq is the minimum beam waist of the incident 
beam located at the entrance pupil of the objective lens. 
The above 1/e- width can be converted to a full width at 
half maximum as 



Axfwhm = 2 v/ln V2 w' ~ 0.375 



NA n 



(16) 



where NA e ff = Wo/ f. 

We can also derive the above result from the I integrals 
for general cases. From Eq. Q, the field distribution in 



2. Longitudinal spot size (Azpwhm) 

In Gaussian optics, the Rayleigh range zq in the region 
of focus is given by 



TTWq' 2 A ( f X 



(20) 



A 7T \UIo, 

The FWHM value in the z direction is just twice of the 
Rayleigh range. 

Azfwhm = 2 (-) (-L) ~ 0.6366—^. (21) 
\irj \w J NA^ ff 

We can also derive Eq. l|2T)l from Eq. Q: 



Jo (r, = 0) 



e o ($) Vcos $ sin 1} ( 1 + cos$) 



x exp(ifcr cos i?) d-d, 
h(r,9 = 0) = Q = I 2 (r,9 = 0). 



(22) 
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Again, the integrand is significant only when i9 < 
Wq/J <C 1, and thus 



Iq oc y exp(— x 2 )a: exp ■|zA:r 

1 Z" 00 

= — exp(ifcr) / exp(— q) exp 

2 Jo 



--kr(w /f) 2 q 



cx 1 + i- 



. /cru> 2 
2F 



r/.r 
dq 
(23) 



and thus the intensity distribution is proportional to a 
Lorentzian 



\Io\ 2 oc 



1 



(2f 2 Y 



(24) 



from which we obtain Azfwhm = 4/ 2 /fcwQ identical to 
the one in Eq. f2l"|) . 



D. Application to NA=0.4, 0.75 and 0.95 

In Fig.^J theoretical Axfwhm and Azfwhm values for 
NA=0.4, 0.75, 0.95, respectively, are plotted. The solid 
lines represent FWHM's calculated from Eqs. I@J and 
©. The dash-dotted lines in Figs. EIa),(c) and (e) are 
calculations based on Eas. p6|l and those in Figs.Q{b),(d) 
and (f) are given by Eq. 1|21|) . Similarly, the dashed lines 
are calculated from Eqs. © and (|T3|) . These dash-dotted 
and dashed lines represent two limiting cases, wq/R <C 1 
and wq/R^$> 1, respectively, of the general curves which 
are valid for any Wq/R values for given NA's. 
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FIG. 4: Dependence of transverse and longitudinal FWHM 
values, Aifwhm and Azfwhm, respectively, on wq of the 
incident Gaussian beam. Vertical dotted lines indicate wo = 
R/2 and R. (a)-(b) NA=0.4, (c)-(d) NA=0.75, and (e)-(f) 
NA=0.95. 



III. EXPERIMENT 
A. Quantity to be measured 

In order to measure the energy flux or the z-component 
of the Poynting vector associated with the field distribu- 
tion near the focal plane, we place a sub-micron pinhole 
at various positions and measure the light power trans- 
mitted by the pinhole. By scanning the pinhole trans- 
versely and longitudinally, we can map out the distri- 
bution of the energy flux. Since the pinhole is made of 
a conductor, the field distribution near the pinhole is 
slightly modified. However, we assume that the effect of 
the interaction between the pinhole and the field on the 
measurement of the energy flux is negligible. The validity 
of this assumption will be discussed in the next section. 

The resolution of an optical microscope is determined 
by the electric field distribution in the focal region of 
the objective lens since samples response to the electric 
field of an illumination light [l6l |. In our experiment, 
however, we measure the energy flux or the z-component 
of the Poynting vector associated with the electric field 



distribution as mentioned above. According to our the- 
oretical investigation, the smallest FWHM of the elec- 
tric field distribution is about the same as that of the 
z-component distribution of the Poynting vector S z in 
the focal region up to NA ~ 1 under our experimental 
conditions. The difference between those two FWHM's 
is about 8%, except for the usual difference, i.e., the de- 
tailed structure in the electric field distribution elongated 
in the incident polarization direction Therefore, the 
transverse and longitudinal FWHM's of S z distribution 
well approximate those of the electric field distribution 
and thus they can be used as measures of the optical reso- 
lution associated with the illumination part of an optical 
scanning microscope. 



B. Experimental Setup 

Our experimental setup is shown in Fig. |S| A He-Ne 
laser (632. 8nm) with x-polarization was first incident on 
a spatial filter, and then expanded and collimated to a 
Gaussian beam with a beam waist wq. Its profile was 
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FIG. 5: Experimental setup for measuring the profile of the 
beam focused by an objective lens. LI, L2, L3: lenses, BS: 
beam splitter, TS1: translation stage driven by a closed- 
loop-feedback stepper motor, TS2: translation stage driven 
by closed-loop-feedback PZT actuators, OL: objective lens, 
CCD: charge-coupled device detector, PI, P2: pinholes, C: 
condenser, PMT: photomultiplier tube, and Al, A2, A3: scan 
control voltage signals from an analog-digital converter board 
on a personal computer. Signal Al controls the z translation 
of the objective lens and signals A2 and A3 control the x, and 
y translation of the pinhole stage. A spatial filter is formed 
by LI, PI, and L2. 



measured by a motorized beam profiler. An objective 
lens was mounted on a xyz-translation stage with its z 
coordinate scanned by a step motor in a closed feedback 
loop (model M-126.PD from PI Ltd.) and thus it could 
be coarse-positioned manually and fine-scanned by the 
step motor with 0.125 /mi resolution in the z direction. 

Infinity-corrected microscope objective lenses with 
NA=0.4, 0.75 (both from NIKON) and 0.95 (from 
OLYMPUS), respectively, were used. The spherical aber- 
ration coefficients of the objective lenses were measured 
with a Twymann-Green interferometer (Zygo) and the 
results are 0.44A, 0.21A, 0.43A for NA=0.4, 0.75, 0.95, 
respectively. In FWHM measurement to be presented 
below the error caused by these values of spherical aber- 
ration is estimated to be negligible, as small as 0.1% or 
less. 

A pinhole (see Fig. with a diameter of (0.50 ± 0.05) 
/tm, which served as an intensity probe, was mounted 
on a translation stage driven by piezoelectric transducer 
(PZT) stacks in a closed feedback loop for scanning in 
the x- and y directions. Typical stroke errors of these 
PZT stacks were less than 0.1% of their stroke ranges. 

The pinhole was made with the electron-beam etching 
technique on a thin Ti:sapphire substrate with a gold 
coating layer of 200 nm thickness. The gold layer in 
a circle of 0.5 micron diameter was removed to form a 
pinhole. 

The light transmitted through the pinhole was de- 
tected by a photomultiplier tube and the signal was digi- 
tized by a data acquisition board in a computer as a func- 
tion of the pinhole position. A resulting image amounted 
to a 200x200 array of pixels. 



FIG. 6: Scanning electron microscopy image of the pinhole 
(4> = 0.5 ± 0.05/tm) used as an intensity probe in our experi- 
ment. 




FIG. 7: Observed zz-profile in the focal region for NA=0.95 
objective lens. The image covers a scan area of 2.5 /im x 6 
/im. 



IV. RESULTS 

We measured the intensity distribution for a Gaussian 
beam with an initial beam waist of u>o=0.57, 0.97, 1.59, 
2.57, 3.1, 3.58, and 5.88 mm. From the measured inten- 
sity distribution in the xz meridional plane, we deter- 
mined FWHM's in the x direction (Asfwhm) and in the 
z direction(A,2FWHM)- 

For instance, the intensity profile created by an objec- 
tive lens with NA=0.95 for an input beam of u>o=5.88 
mm is shown in Fig. \7\ Since the entrance pupil radius 
R of the objective lens was 1.71 mm, we can consider 
the incident beam as a plane wave. The x-z profile cor- 
responded to an actual area of 2.5 /mi x 6 /mi. The 
measured x- and z-FWHM values were 0.4 /mi and 1.03 
/im, respectively. 

The time-averaged z component of the Poynting vec- 
tor in the near focus was calculated from Eq. (JSJ. To 
compare experiment with theory, we assumed that the 
total amount of light detected by the PMT through the 
pinhole was proportional to the convolution of the z com- 
ponent of the Poynting vector with the pinhole opening. 
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S z (x,y) = / S z (x',y')P(x-x',y-y')dx' dy 1 (25) 



where P(x, y) is an aperture function for the pinhole. 
This assumption is equivalent to saying that the possi- 
ble field distortion by the conducting surface of the pin- 
hole substrate does not affect the amount of energy flow 
through the pinhole much so that we just integrate the 
surface-normal component of the unperturbed Poynting 
vector calculated for the absence of the pinhole over the 
aperture function of the pinhole. 

The dependence of the measured x- and z-FWHM val- 
ues on the input Gaussian waist wq for NA=0.4, 0.75 
and 0.95 are shown in the Fig. [5] where (a) and (b) are 
for NA=0.4, (c) and (d) for NA=0.75 and (e) and (f) 
for NA=0.95, respectively. FWHM values obtained from 
Poynting vector S z (x, y) and convoluted Poynting vector 
S z (x, y) are represented by solid and dashed lines, respec- 
tively and experimental results are represented by square 
dots. The spherical aberration of the object lenses was 
included in the calculation of S z . The agreement between 
experiment and theory is reasonably good. 

The smallest (rFWHM, z-FWHM) values measured 
in the experiment are (0.86 /im, 7.62 /mi) for NA=0.4, 
(0.48 fjm, 1.79 /xm) for NA=0.75 and (0.40 /im, 1.03 fim) 
for NA=0.95. These values are in good agreement with 
the convoluted FWHM values except for the x-FWHM 
values for NA=0.75 and 0.95 (see Figs. Etc) and (e)), 
for which the observed FWHM is slightly smaller than 
the convoluted FWHM but larger than the un-convoluted 
FWHM. 



DISCUSSIONS 



Effect of the interaction between the pinhole 
and the electric field 



The experimental results summarized in Fig. [S] show 
that the observed FWHM is smaller than the FWHM of 
the pinhole-convoluted S z distribution when the focused 
beam spot size is comparable to or smaller than the pin- 
hole size. For such small focal beam spots, the pinhole 
seems to behave as a smaller pinhole for light transmis- 
sion. This phenomenon appears to be caused by the in- 
teraction of the electric field and the conducting surface 
of the pinhole. The distortion of the field distribution 
near a conducting structure like a pinhole is usually in 
the sub-wavelength scale and thus it can be neglected if 
the range of the field distribution is much larger than 
the wavelength. If the range of the field distribution is 
in the sub-wavelength scale, as in the case of x-FWHM 
for NA=0.75 and 0.95, the field distortion effect could be 
non-negligible. Our numerical simulation supports this 
reasoning. Nonetheless, the effect is still small, amount- 
ing to at most 20% with respect to the simple-minded 
convoluted FWHM, under our experimental conditions 
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FIG. 8: Dependence of transverse (x) and longitudinal (z) 
FWHM values on wo of the incident Gaussian beam. Uncon- 
voluted FWHM's obtained from Eq. © are represented by 
solid lines whereas the convoluted FWHM's given by Eq. ()25fl 
are drawn as dashed lines. Experimental results are marked 
by square dots with error bars. Independently measured 
spherical aberrations were included in the calculations. Ver- 
tical dotted lines indicate wo — R/2 and R. (a)-(b): NA=0.4, 
(c)-(d): NA=0.75, and (e)-(f): NA=0.95. 



and thus our previous assumption of neglecting this effect 
could be well justified in the first-order approximation. 



B. Range of input beam waist for acceptable focal 
spot size 

When wq is equal to the radius of the entrance pupil 
R, the resulting theoretical x-FWHM values are larger by 
about 10% for all three NA's than the ultimate FWHM 
values, which occur when wo ^S> R- For NA=0.4, 0.75 
and 0.95, i?=4.0 mm, 3.0 mm and 1.71 mm, respectively. 
The difference between the theoretical z-FWHM's and 
the ultimate z-FWHM's are 7.3%, 4.8% and 3.5% for 
NA=0.4, 0.75 and 0.95, respectively, in this case. When 
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wo = -R/2, the difference increases to about 50% for x- 
FWHM for all three NA's and to 88%, 65% and 54% 
for NA=0.4, 0.75 and 0.95, respectively, for z-FWHM. 
In other words, the FWHM of the Gaussian beam at the 
focus docs not increase much until the input beam waist 
is reduced below the half of the entrance pupil radius of 
the objective lens. Our experimental results support this 
observation. 



C. Effect of quasi-Gaussian input beam 

There exists a small discrepancy between theory and 
experiment in the regime of wq < R/2. It is attributed to 
the use of an imperfect Gaussian beam as an input beam. 
In real experiments, the laser beam is not a perfect Gaus- 
sian beam. In order to make it close to a true Gaussian 
beam, spatial filtering of the laser beam is performed. For 
most of laser applications, a single-pass spatial filtering 
is more than enough. We have found, however, in our ex- 
periment with wo <C R a single-pass spatial filtering was 
far from adequate since in this case the entire beam pro- 
file of the input beam determines the field distribution in 
the focal region. 

Although a spatial filtering process can eliminate most 
of the asymmetric structures in the spatial distribution of 
an incident beam, the resulting mode distribution tends 
to contain small side lobes of Airy disk type. Unless these 
side lobes are thoroughly eliminated by a succession of 
extensive spatial filtering, the resulting field distribution 
in the focal region becomes significantly broadened and 
distorted from the expected distribution of the true Gaus- 
sian beam. 

According to our numerical simulations, the error in 
FWHM measurement induced by the imperfect Gaussian 
beam may amount to 30% or more and the error is par- 
ticularly considerable in the regime of wo -C R. In our 
experiment, we have used three successive stages of spa- 
tial filtering in order to minimize any deviation from the 
true Gaussian beam. When only a single stage of spa- 
tial filtering was used, we observed about 30% increase 
in FWHM's in most cases. 



VI. SUMMARY AND CONCLUSIONS 

We studied the intensity distribution in the region of 
focus when a linearly-polarized well collimated Gaussian 
beam with a waist of Wo was incident on a high-NA ob- 
jective lens with an entrance pupil radius of R. We first 



theoretically examined the transverse spot size Axfwhm 
and the longitudinal spot size Axfwhm near the focal 
plane for an arbitrary input beam waist wq. We used 
the vectorial diffraction theory of Richards and Wolf 
and calculated a time-averaged Poynting vector in the 
near focus. We then derived approximate expressions for 
FWHM's for two limiting cases, wo -C R and w 3> R, 
and for the latter the approximate expression is in the 
form of Fraunhofer diffraction result although the result 
is obtained for the non-paraxial case. 

In experiments, we varied the initial wo for a given 
NA's of 0.4, 0.75 and 0.95 and measured Axfwhm and 
A^fwhm values. They were obtained by scanning a pin- 
hole of 0.5 /im diameter across the focused beam and by 
measuring the total transmitted light through the pin- 
hole. The results obtained by convoluting the calculated 
Poynting vector with the pinhole were well matched with 
the measured intensity distributions. The smallest mea- 
sured x- and z-FWHM values were 0.40 ^m and 1.03 /jm, 
respectively, for NA=0.95 with A= 632.8 nm. 

For high NA's with resulting focal beam waists com- 
parable to or smaller than A, observed i-FWHM's were 
smaller than those of the energy flux distribution convo- 
luted with the pinhole. This discrepancy is attributed to 
the distortion of the electric field near the conducting sur- 
face of the pinhole. In addition, we observed a small dis- 
crepancy between theory and experiment for wo < R/2, 
which is caused by slight deviation of the spatial distribu- 
tion of the incident beam from that of the true Gaussian 
beam. Extensive multi-stage spatial filtering was used to 
minimize this deviation. 

Finally, we found both theoretically and experimen- 
tally that the FWHM of the Gaussian beam at the focus 
does not increase much until the input beam waist is 
reduced below the half of the entrance pupil radius of 
the objective lens. This result can be used as a practi- 
cal design guideline for scanning microscopy employing a 
Gaussian laser beam as a probe. It has been also noted 
that the spatial distribution of the incident beam has to 
be as close to that of a true Gaussian beam as possi- 
ble via extensive spatial filtering, in order to achieve the 
smallest focal beam spot, particularly when a beam with 
wo < R/2 is used as a scanning probe. 
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